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Abstract 

We study FRW homogeneous cosmological solutions in the bigravity 
extension of the recently found ghost-free massive gravity. When the 
additional extra metric, needed to generate the mass term, is taken as 
nondynamical and flat, no homogeneous flat FRW cosmology exists. 
We show that, when the additional metric is a dynamical field a per- 
fectly accetable FRW solution exists. Solutions fall in two branches. 
In the first branch the massive deformation is equivalent to an effective 
cosmological constant determined by the graviton mass. The second 
branch is quite rich: we have FRW cosmology in the presence of a 
"gravitational" fluid. The control parameter £ is the ratio of the two 
conformal factors. When £ is small, the evolution is similar to GR 
and interestingly the universe flows at late time towards an attractor 
represented by a dS phase. 



1 Introduction 

Recently, there has been a renewed interest in the search of a modified theory 
of gravity at large distances through a massive deformation of GR (see P 
for a recent review) . A great deal of effort was devoted to extend the seminal 
work of Fierz and Pauli (FP) [3j at nonlinear level [2]. The FP theory, de- 
fined at linearized level, is plagued by a number of diseases. In particular, the 



modification of the Newtonian potential is not continuous when the mass m 2 
vanishes, giving a large correction (25%) to the light deflection from the sun 
that is experimentally excluded [4] . A possible way to circumvent the physi- 
cal consequences of the discontinuity was proposed in [5]; the idea is that the 
linearized approximation breaks down near a massive object like the sun and 
an improved perturbative expansion must be used that leads to a continu- 
ous zero mass limit. In addition, FP is problematic as an effective theory. 
Regarding FP as a gauge theory where the gauge symmetry is broken by a 
explicit mass term m, one would expect a cutoff A2 ~ mg~ l = (mMpi) 1 / 2 , 
however the real cutoff is A5 = (m 4 M p ;) 1 / 5 or A3 = (m 2 M p ;) 1//3 , much lower 
than A2 [6]. A would-be Goldstone mode is responsible for the extreme UV 
sensitivity of the FP theory, that becomes totally unreliable in the absence 
of proper UV completion. Recently it was shown that there exists a non 
linear completion of the FP theory that is free of ghosts up to the fourth 
order [7] and avoids the presence of the Boulware-Deser instability [8]. Then 
the propagation of only five degrees of freedom and the absence of instabili- 
ties was extended to the non perturbative level in [9]; this was shown also in 
the Stiickelberg language in [10]. The bigravity extension of this thoery was 
shown to be ghost-free 

Quite naturally massive gravity leads to bigravity. Indeed, any mas- 
sive deformation, obtained by adding to the Einstein-Hilbert action a non- 
derivative self-coupling of the metric g, requires the introduction of an ad- 
ditional metric g. This auxiliary metric may be a fixed external field, or be 
a dynamical one. When g is non dynamical we are dealing with aether-like 
theories; on the other hand if it is dynamical we enter in the realm of bi- 
gravity [12], originally introduced by Isham, Salam and Strathdee [H]. The 
need for a second dynamical metric also follows from rather general grounds. 
Indeed, it was shown in [15] that in the case of non singular static spherically 
symmetric geometry with the additional property that the two metrics are 
diagonal in the same coordinate patch, a Killing horizon for g must also be 
a Killing horizon for g. Thus, it seems that in order that the Vainshtein 
mechanism is effective and GR is recovered in the near horizon region of 
a black hole, g has to be dynamical [16]. In this paper we also show that 
cosmology calls for the bigravity formulation of massive gravity. While in 
the Stiickelberg formulation there is no homogeneous flat FRW solution [T7] . 
(see also for a related work [TSJ Q3JJ |2UJ [2T| ) in the present work we show that 
flat FRW homogeneous solutions do exist in the bigravity formulation. 

In section [2] the formulation of massive gravity as bigravity is reviewed. 
The cosmological ansatz is introduced in section [31 where the structure of the 
modified Einstein equations and the consequences of Bianchi identities are 
studied. Cosmological evolution falls in two branches described in section 0] 
and section [5J In section El the results of the previous sections are extended 
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to the case of spatially curved geometries. Section [7] contains our conclusions. 
The full set of Einstein equations can be found in the appendix. 

2 Massive Gravity and Bigravity 

Any modification of GR that turns a massless graviton into a massive one 
calls for additional DoF (degree of freedom). An elegant way to provide them 
is to work with the extra tensor g^. When coupled to the standard metric 
g^,, it allows to build non-trivial diff- invariant operators that lead to mass 
terms, when expanded around a background. Consider the action [12] 

S = J d A x^gnM 2 pl 1l+^ [M 2 pl (U - 2m 2 V) + L matt ] , (1) 

where R(gi) are the corresponding Ricci scalars and the interaction potential 
V is a scalar function of the tensor X£ = g^ a g a v Matter is minimally coupled 
to g and it is described by L matt . The constant k controls the relative size 
of the strength of gravitational interactions in the two sectors, while m sets 
the scale of the graviton mass. The action ([T|) brings us into the realm of 
bigravity theories, whose study started in the '60 [13] . An action of the form 
(00) can be also viewed as the effective theories for the low lying Kaluza-Klein 
modes in brane world models [12] ■ The massive deformation is encoded in the 
non derivative coupling between g^ and the extra tensor field g^. Clearly 
the action is invariant under diffeomorphisms, which transform the two fields 
in the same way (diagonal diffs). Taking the limit k — >■ oo, the second metric 
decouples, and gets effectively frozen to a fixed background value so that the 
"relative" diffeomorphisms are effectively broken, as far as the first metric is 
concerned. Depending on the background value of g yw one can explore both 
the Lorentz-invariant (LI) and the Lorentz-breaking (LB) phases of massive 
gravity [13]. When the second metric is dynamical this is determined by its 
asymptotic properties, as discussed below. In this case notice that g^ u is 
determined by its equations of motion (for any finite M p j) so that we will be 
working always with consistent and dynamically determined backgrounds. 
The role played by g^ u is very similar to the Higgs field, its dynamical part 
restores gauge invariance and its background value determines the realization 
of the residual symmetries. 

The modified Einstein equations can be written as3 

E " + Q^=2W i T " ( 2 ) 
kE£ + Q£ = 0- (3) 

1 When not specified, indices of tensors related with g(g) are raised/lowered with g(g) 
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where we have defined Q± and Q 2 as effective energy-momentum tensors 
induced by the interaction term. The only invariant tensor that can be 
written without derivatives out of g and g is X£ = gi a g av [12] . The ghost free 
potential [71 ITT]R V is a particular 4-parameter scalar function of Y£ = (y/~X)£ 
given by 

4 

V = J2a n V n , n = 0...4 (4) 

n=0 

where The V n are the symmetric polynomials of Y 

V = l V 1 =t 1 , V 2 = t?-t 2 , ^ 3 = r 1 3 -3r 1 r 2 + 2r 3 , 

(5) 

V4 = r\ — 6 r 2 r 2 + 8 T\ t 3 + 3 r 2 — 6 r 4 , 

with r„ = tr(y n ). In [11] it was shown that in the bimetric formulation the 
potential V is ghost free as in the Stiickelberg formulation. We have that 

QiiJ = m 2 [ V ^ - (V y)£] (6) 

Q 2 M = m 2 q -l/2 {V > Y) ^ (7) 

where (V'Y U = dV/dY» and g = detX = det(g) / det(g). 



3 Ansatz, Equations of Motions and Conservation Laws 

For simplicity, here we consider the case of spatially flat geometries, the non 
flat case is discussed in section |6j 

In general we cannot set both metrics in diagonal form and preserve 
homogeneity. For instance, take the following form for g and g: 



(8) 



ds 2 = a 2 (t) (-dt 2 + dr 2 + r 2 dQ 2 ) 
ds 2 = u 2 {t) [-c 2 {t) dt 2 + 2D(t) dt dr + dr 2 + r 2 dQ 2 ] 

It is convenient to define 

( = ^. ( 9 ) 
a 

If D 7^ 0, the metric g is still homogeneous and the hypersurface of homo- 
geneity is t — const; on the other hand the metric g is homogeneous with 
respect to a different time slicing. In addition, Qi, Q 2 are not diagonal and 
Qi/2t Qi/2g — as a result Qi and Q 2 are not energy momentum 

tensors (EMT) for a homogenous perfect fluid. In addition, using Bianchi 
identities and the equations of motion one can show that there is no physical 



2 A very similar potential having the same form but with X instead of X 1 / 2 was con- 
sidered in l22l. 
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solution when D(t) ^ 0. Indeed, from the fact that the Einstein tensor for g 
is diagonal, also Qi should be diagonal and with D^O one finds that £ = £ 
is constant with 6a 3 £ 2 + 4a 2 £ + a\ — Furthmore, Bianchi identities can 
be fulfilled only if a 2 + 3 £ a 3 = 0. Finally, in order to solve the tt component 
of Einstein equations for the metric g for any value of r, we are forced to set 
D — 0. Thus, we only need to consider the case D — 0. 

The potential part of the action is separately invariant under diff and 
thus it gives two set of Bianchi identities V v Qi^ = V 1 = which are 
equivalent to the following conservation 

±(a>P,)+p4/ = 0, (10) 

where 



Pg = £^ (3 a 3 £ 3 + 3 a 2 £ 2 + | a x £ + y 



m 



2 



3a 3 c£ 3 + a 2 (2 C +l)£ 2 + ^( C + 2)£ 1 



47rG L * v ^ 2 v ^ 2 J 
Eq. ffTUl) is equivalent to the following condition 

m 2 (6a 3 £ 2 + 4a 2 £ + a!) (cua'-au 1 ) = 0. (12) 

The conservation of the matter EMT gives 

j {a 3 p m ) +Pmj f a 3 = 0. (13) 

In order to determine the three functions of time a, c, a; we need a set of 
three independent equations. As in GR, one can show that for both metrics 
the time-time plus the space-space component of the equations of motion 
leads to the Bianchi identities. As a result, also for massive gravity we can 
take the conservation equations (U2H13p . together with 

H 2 ( Pm + Pg ) , (14) 

H 2 ?7l 2 

: (ax C 3 + 6 a 2 C 2 + 18 a 3 C 1 + 24 a 4 ) , (15) 



c 



2 3k 



as independent equations. We defined the Hubble parameters in conformal 
time H = a' /a 2 and H u = u' /u 2 . The full set of equations of motion is given 
in Appendix lAl 



3 Notice that the constraint coincides with the one found for branch one solutions below, 
see sect. |U 

4 One can easily show that diff invariance also implies that V 1 V Q\ U = <==>■ V ' V Q2 U = 0. 
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Starting from eq. (|T2|) the cosmological evolution can be classified in two 
branches: either there is an algebraic constraint for £ 

6a 3 £ 2 + 4a 2 £ + ai = 0, (16) 

or 

- = c~. (17) 

uj a 

Notice that, by definition c(t) is the positive root of c 2 and thus u' and a' 
must have the same sign. Since in our universe a' > 0, we will consider only 
the case oj' > 0. 



4 Cosmology: Branch one 

Let us consider the case when the conservation equation (fT2l is solved as 
in ffTB"]) . so that clearly £(t) = £ is a constant satisfying the condition 

6a 3 f + 4a 2 £ + ai = (18) 

and E w = iff- 1 . 

From (fill) we have 

H 2 = ^ Pm + ^ [oo - 6f (2a 3 f + a 2 )] . (19) 

As a result, in this branch the effect of the massive deformation is to induce 
an effective cosmological constant 

Ae// = ^ " ^ ( 2fl3 ^+ fl2 )] • ( 2 °) 
For what concerns the second metric, (]15[) turns into an algebraic equation 
for c(t) : 

c 2 = 3kH2 (21) 

2m 2 [6£(2a 4 £ + a 3 ) + a 2 ] ' 

Summarizing, branch one is equivalent to the usual FRW cosmology with a 
cosmological constant. This is not very surprising. Let us remind that in the 
bimetric formulation, when one expands around flat space at the quadratic 
level, there are two spin-2 modes: one is always massless by diff-invariance, 
the other has Pauli-Fierz mass, see for instance [121 03] • The point here is 
that condition ( TT8]) leads to a zero quadratic Pauli-Fierz mass. The condition 
is the same found in [16] in the branch where no gravity modification are 
present. Notice that, even if at the quadratic level the potential has no 
effect, reappears in the theory at higher orders in perturbation theory. 

Finally, let us comment on the k — > 00 limit which corresponds to freezing 
the second metric. From eq. (|15[) one gets that H u = 0, but in the branch we 
are dealing with, £ is constant and so also H = 0. As result, for the branch 
one no flat FRW cosmology exists when the second metric is non dynamical. 
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5 Cosmology: Branch two 

The Bianchi identities ( 1T7|) in such a branch give 



c=§^- (22) 
Inserting the expression of c in (IT51) we have 

m 2 [6£(4a 4 £ 2 + 3a 3 £ + a 2 ) + aj -3k£H 2 = 0. (23) 
If we solve this equation for if 2 and compare with ( Tl4|) . we find 



m 2 



f 2 a 2 + H ai + - — - + 2 a 3 £ - — 

\ K / \ AC J 6 K t, K 6 

= p m ^ 

For matter with an equation of state p m = wp m , we have that p m = p a~ 3( - 1+w \ 
Thus, eq. ( |24p can be solved for £ in terms of the conformal scale factor a(t). 
The problem of the cosmological evolution is reduced to FRW cosmology 
in the presence of matter plus an additional form of "gravitational" mat- 
ter. Though the solution of f l24|) can be written in a closed form, it is not 
particularly illuminating. 

Alternatively, from the definition of £ and from (|22p we have 

e=(c-l)aH£, (25) 

and from (j2l|) we can express a in terms of £: a = f a (C,,w). Finally, from 
(158]) and fl60|) (see appendix [A} c can be expressed in terms of £ 

c = 1 + (1 + ^) / c (0 , (26) 

fdO is a function whose explicit form is not needed. At this point, the 
velocity field of £ is function of £ only 

e = a + w) m, «o fait, w) e , (27) 

and it has two fixed points 

= c = lwith/ c (£) = p m ^0 
£ = £ , with tt>) = generic p,; 



(28) 



'in 



While the fixed point with zero Hubble parameter is not cosmologically very 
significant, the c = 1 case instead corresponds to a De Sitter phase. 
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Far from the fixed points the main features can be deduced by a qualita- 
tive analysis of the solutions of (124)) . The variable £ is dimensionless and natu- 
rally is a function of G p m j m 2 . In the early time, when G p m j m 2 > 1 we have 
two possible regimes: small £ oc 0{m 2 /G p m ) and large £ oc 0{G p m /m 2 ). 



Large £ 



The regime of large £ is determined by the leading positive power of £ present 
in eq. ( |24l) . For simplicity, here we report only the leading results, the cor- 
rections to them are O (m 2 /G p m ). For a 3 7^ the Leading term is £ 3 and a 
solution exists only if a 3 < 0, with 

'8irG p r , x 1/3 



Pm + Pg 
W eff 



&\a 3 \m? 



4 G4 — K 0,2 
3 d 3 K 



k V ir G a 2 



Pg+Pn 
Pg + Pr, 

—W . 



2w-l 



(29) 

(30) 

(31) 
(32) 



When a 3 = 0, the leading term is £ 2 and a solution exists if 4 a 4 — k a 2 > 0, 
with 

1 l/2 

; (33) 



8ttGk Pr, 



m? 6(4 (24 — Ct2 ft) 

4a 4 



(4 04 - k a 2 
3w + 1 



-Pm ; (34) 

(35) 
(36) 

Finally, if a 3 = and fta 2 - 4a 4 = 0, the leading term is £ and a large £ 
solutions exists when a\ < 0, with 

8ttG 



Pm + Pg 

W eff 
C 



3 ai m 2 



Pm i 



16 7r G a 2 2 
3 a\ m z 

2w + l; 
-(3w + 2). 



(37) 

(38) 

(39) 
(40) 
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The large £ regime is characterized by non standard cosmology and/or nega- 
tive values of c. Except for the case where the quadratic term is dominating, 
the effective equation of state w e ff largely deviates from w, but even in this 
case c < when w > 0. As a result, we conclude that the large £ regime is 
not a physical one. 



Small £ 

Solutions with small £ are present only when a\ > and in this case 

t = V~F^ ' 41 

8 7T G K p m 

Pm + Pg = Pm] (42) 
W eff = W, (43) 

c = A + 3w. (44) 

For simplicity, we give only the leading terms, the corrections are 0(m 2 /G p m ) 
and can be sistematically computed. In this regime cosmology is standard: 
once the matter is so diluted that p m is negligible in (124"]) . the system falls in 
the fixed point region and £ is almost constant; the universe enters in a late 
time dS phase. 

One can solve perturbatively ( fl4l) to find a. For a radiation dominated 
universe we fined 

t a m 2 t 5 mH* (alK + 20aj) , / 3 \ 1/2 

ait) ~ T + ^0^ + 1080^ + ' ' ' ' t0 ~ ) ■ (45) 

In the case of matter dominated universe 



2 „ ^,2+8 ^4 + 14//i^2^ | /in^2N / o \ V 2 



t 2 a m 2 t 8 m 4 t 14 (4a 2 1 fi; + 49a 2 , 
a(t) = + — — =- + i — 2 77: ^ + 



t 2 84 tg 30576 V 27r Gpo 



(46) 



For instance, taking po of order of the critical density today and m ~ 10 _33 e^/ 
we have that during the radion era £ ~ 10 (1 + z) ~ 4 << 1. 

Let us also discuss the k — > 00 limit in this branch. When the second 
metric is non dynamical, from eq. f)23p we get that H — > 0. Thus, also in this 
branch no flat FRW cosmology exists with just a single dynamical metric. 



5 As an example, in the following expressions we give the leading and the next to leading 
correction to a. 
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6 FRW with Spatial Curvature 

Let us discuss the case with spatial curvature (see [20J for the case with frozen 
second metric). The metrics take the form 



ds 2 = a 2 (t) (-dt 2 + - Ts + r 2 dtf) 

\ (1-^r 2 ) J 



7 2 



ds =uj 2 (t) 



dr 2 < 47 > 



-c 2 (t)dt 2 + - -— + r 2 dn 2 

(1 - k 2 r l ) 



In the presence of curvature the conservation law for Q\ and Q 2 take a 
different form 

(F-l) [2c£ (3a 3 i + a 2 )+2a 2 £ + a 1 ] = 0, (48) 

F(r) £ {[2 £ (9 a 3 £ + 4 a 2 ) + aj H u - 4 c H (3a 3 £ + a 2 ) } - 

cH (6a 3 £ 2 + 8a 2 £ + 3ai) + 2f (2a 2 f + ai) = 0, (49) 



r 2 1 \ 1/2 
ki r z — 1 x 



where 

F M = (S-lJ ■ (50) 

When k\ 7^ fc 2 , F is different from one and the Bianchi identities must hold 
for any value of r; this gives three nontrivial relations. From (l4*8|) . we can 
solve for c 

2a 2 £ + di 



6a 3 £ 2 + 2a 2 £' 
then inserting c in ( T4*9]) we get that 



(51) 



2ff(2a 2 £ + a 1 ) 

if (2a 2 e + a 1 )(6a 3 e 2 + 4a 2 e + ai) 2 = 0. (53) 

The solution 2 a 2 £ + ai = would lead to c = and to a degenerate second 
metric. Thus, either H = H w = or £ = £ = const, with 

6a 3 £ 2 + 4a 2 £ + ai =0. (54) 

In the special case a 2 = a 3 = we get the nonphysical solution H = H w = 0. 
Thus, the only interesting solutions of Bianchi identities, when k\ ^ k 2 , are 
the ones with £ = £ satisfying eq fl54l) . Notice that for £ = £ we have that 
c—1. The further constraint coming from the equations of motion gives 

k 2 - h 8tiG , rr , 
— + -z~p m + A = , (55) 

(1 o 
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where A is a function of £. The previous equation is algebraic in a and it gives 
a =constant and H = 0. Thus, in order to obtain a reasonable cosmology we 
are forced to set k\ — — k c and we are back to the Bianchi identity (fT2l) . 
The analysis of the previous sections is the same, except for the presence of 
the curvature term. 

Let us now discuss the differences with the Stiickelberg approach where 
g is non dynamical. Formally, the non-dynamical limit would correspond to 
hi — y oo. To make contact with the existing literature, we take g equivalent 
to the Minkowski flat metric implying clearly that c and u cannot be arbi- 
trary. Imposing that the Riemann curvature tensor of g vanishes we find the 
conditions 

fc c <0, c=^=. (56) 

Therefore, FRW cosmology with frozen second metric exists only with a 
negative non- vanishing spatial curvature [17]. When ( )56|) holds, Bianchi 
identities can be realized only within branch one, leading to eq. (fT8j) . 

7 Conclusions 

In this paper we studied FRW cosmological solutions of the bigravity exten- 
sion of the recently found ghost free massive gravity theory. It was shown 
in ref. [TTJ 120] that when the auxiliary metric needed for the formulation of 
the theory is taken to be non dynamical, no flat FRW cosmology exists. On 
the other hand, as we showed here, flat FRW solutions are allowed when the 
second metric is promoted to a dynamical field. This result, together with 
the analysis of spherically symmetric solutions [16], indicate that bigravity 
is more than just a tool; it is an important ingredient in the formulations of 
a physically acceptable theory of massive gravity. The cosmological evolu- 
tion is governed by Bianchi identities (BI). Depending on the way they are 
realized, cosmology falls in two separate branches. When the BI are imple- 
mented by an algebraic constraint on the ratio £ of the two scale factors, we 
find standard FRW cosmology, in the presence of a cosmological constant 
proportional to the graviton mass scale. When BI are implemented by a 
differential relation between the Hubble parameters of the two metrics, cos- 
mology is instead richer. Cosmological evolution at the early times, where 
Gp m /m 2 3> 1, is controlled by the parameter £. In the presence of matter 
with an equation of state w, the large £ regime leads to unphysical solutions 
with c < 0. On the contrary, the small £ regime is safe and reproduces the 
standard early FRW cosmology; at late time the universe sets in a dS attrac- 
tor region. The model is interesting and though there are a number of free 
parameters in the potential, it is surprisingly predictive. The next step is to 
study perturbations to check whether massive gravity cosmology is viable. 
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We will report these results in a future work. 



Note added 

During the completion of this work, a study on the same topic appeared [23 
Our results seem to agree with [23], except for some exotic case. Almost at 
the same time of our paper, another paper [2^ on the same subject was 
announced in the arXiv. 



A Modified Einstein Equations 

For completeness we collect the full set of equations and we restore the pres- 
ence of the spatial curvature k. The non vanishing equations for the metric 
g are the tt component and the space-space components and read 

3k 

m 2 [6a 3 £ 3 + 6a 2 £ 2 + 3a^ + a l + 8vrGp m = — + 3H 2 ; (57) 

cr 

m 2 [6a 3 c£ 3 + 2a 2 (2c + 1)£ 2 + a x (c + 2)f + a ] - 8?rGK 

= - 2 + — + 3H 2 . (58) 
a 2 a 

For the metric g the structure is the same and we have 

, ( a\ 6a 2 18<23 ni \ 3k k 3kH 2 



2 / ai 4a 2 6a 3 2a 2 12a 3 

+ ^ + (60 ) 

Except for section El in order to de-clutter formulae we have not included 
k; as we can see from the previous equations the curvature can be easily 
restored without effecting the analysis. A linear combination of eq. ( 157]) and 
its time derivative together with (|58|) gives ([12]) and the conservation law for 
matter. The same is true for fl59l) and fl60l). 
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